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NON-EXISTENCE  OF  SNAP- THROUGH  FOR  CLAMPED, 
SHALLOW  ELASTIC  ARCHES  SUBJECTED  TO  IMPULSIVE  LOADING1 


B.  Vahidi2 


ABSTRACT 

The  problem  of  dynamic  snap- through  instability  of  the  elastic 
clamped  shallow  arches  has  been  re-examined.  The  governing  equation 
for  the  corresponding  static  problem  is  solved,  and  it  is  proved  that  for 
ideal  impulsive  loading  there  exist  nc  other  symmetric  or  assymmetric 
static  equilibrium  configurations  aside  from  the  original  undeformed 
position.  The  usual  definition  of  the  snap- through  adopted  i„  this  paper, 
requires  existence  of  more  than  one  static  equilibrium  configuration. 
Hence,  it  is  concluded  that,  contrary  to  previous  analyses,  clamped 
shallow  arches  under  impulsive  loading  do  not  snap- through. 
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NOMENCLATURE 


2a  =  span  of  arch 
A  =  cross-sectional  area  of  arch 
A1.A2.A3.A4  =  constants 
E  =  Young's  modulus 

H  =  rise  of  arch 
I  =  moment  of  inertia  of  arch 

p  =  mass  of  arch  per  unit  length 
t  =  time 

v/q  =  original  shape  of  the  arch  measured  from  baseline 
w  =  deformed  shape  of  the  arch  measured  from  baseline 
r  =  horizontal  coordinate 
q  =  external  load 

x  =  —  ~  nondimensional  horizontal  coordinate 

a 


1 

y0  "  2 


nondimensional  original  arch  shape 


=  nondimensional  deformed  arch  shape 
rise  parameter 

=  nondimensional  load 


nondimensional  time 


1.  INTRODUCTION 


The  problem  of  the  dynamic  snap- through  stability  of  shallow 
elastic  aiches  has  been  studied  by  several  investigators  [1,  3,  5]  . 

These  studies  have  been  mainly  c'nfined  to  the  arches  subjected  to: 

(1)  step-loading  where  a  constant  load  is  suddenly  applied  at  zero  time 
and  maintained  thereafter;  (2)  ideal  impulsive  loading,  v’hich  is  assumed 
to  be  equivalent  to  imparting  an  initial  velocity  to  the  arch. 

The  snap- through  stability  of  elastic  hinged  arches  subject  to 
impulsive  loading  has  been  considered  in  a  comprehensive  manner  by 
Ksu  [5,  6]  .  He  has  derived  expressions  for  the  static  equilibrium  con¬ 
figurations  and  has  obtained  sufficient  conditions  for  stability.  For  the 
dynamic  analysis  of  the  clamped  arches,  however,  only  approximate 
solutions  have  been  developed.  Simitses  [13  assumes  a  deflection  shape 
in  the  form  of  a  series  which  includes  two  symmetric  and  one  antisym¬ 
metric  mode.  Then, by  minimizing  the  expression  for  the  Hamiltonian, 
he  obtains  a  set  of  ordinary  differential  equations  from  which,  through 
use  of  the  potential  surface  method,  he  obtains  upper  and  lower  bounds 
on  the  critical  impulse.  Humphreys  [2]  also  has  assumed  a  deflection 
shape  in  the  form  of  a  series;  however,  he  retains  the  first  six  terms  ol 
the  series  expansion.  The  solution  is  then  obtained  with  the  aid  of  an 
analog  computer.  In  [2]  there  are  also  presented  experimental  data  in 
an  attempt  to  verify  the  numerical  results. 


Numbers  in  square  brackets  designate  references  listed  at  the  end 
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In  the  present  paper,  the  snap-through  of  circular  arches  with 
clamped  ends  subjected  to  ideal  impulsive  loading  is  re-examined. 

Throughout  the  analysis  it  is  assumed  that  the  arch  is  perfect,  shallow  and 
elastic.  To  obtain  possible  static  equilibrium  configurations,  the  cor¬ 
responding  static  equations  under  zero  load  are  solved  exactly,  in  a  manner 
suggested  by  Schreyer  and  Masur  [4]  .  It  is  then  shown  that  these  equations 
have  no  solutions  for  any  real  values  of  axial  load  other  than  zero,  indicating 
that  the  only  static  equilibrium  configuration  for  the  arch  is  the  original 
undeformed  position.  It  follows  from  the  accepted  definition  of  snap- through 
that  snap- through  cannot  occur  if  there  is  only  one  equilibrium  configuration. 
[1,  3,  5,  6]  .  Hence,  it  is  concluded  that  the  shallow  circular  arch  with  clamped 
ends  under  impulsive  loading  does  not  snap- through. 
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2.  ANALYSIS 


In  the  notation  used  in  [3]  ,  and  with  reference  to  Fig.  1,  the 
equations  of  motion  for  an  elastic  shallow  arch  can  be  written  as: 

a4(w-wrt) 


EI - 2.'  +  N  &S  4  p 

cr4  Sr2  at2 


+  q  =  0 


where  N  is  the  axial  force  in  the  arch  given  by  the  relation, 

a 


N  = 


AE 

4a 


im  -  ter] 

-  a 


dr 


The  above  equations  can  be  transformed  to  nondimensional 
form  by  the  introduction  of  the  following  dimensionless  quantities: 

I  U/2 


wo  s2(t) 


w 


I  \l/2 


r 

x  s  — 

a 


/2 


’■HW 


.  ^  /AjI/2 

^  O  T7T 


8“  = 


2EI  \  I  / 

a2  N 
EI 


Equations  (1)  and  (2)  in  dimensionless  form  are  then 

54(y-y0) 


+ 


«■  P‘ 


+  Q  -  0 


Bx4 


5x‘ 


(1) 


(2) 


and 


B‘  =  ^  [(^r)  -  (S)2  ] d* 

Boundary  conditions  fcr  an  arch  clamped  at  both  ends  are 
given  below 


w  (-a,  t)  =  0 

(5a) 

w  (a,  t)  =  0 

(5b) 

d  w  .  t  _ 

Sr  '1  ~ 

dw 

-37  <-a> 

(5c) 

—  (a  t'  = 

dw 

dr  <«> 

(5d) 

Now,  let  the  initial  shape  of  the  arch  be  given  by 


or,  in  a  dimensionless  form, 

yQ  ~  X  (1  -  xs)  (7) 

The  boundary  condition  (5a-d)  in  dimensionless  form,  with 
yQ  ,  given  by  (7)  can  now  be  written  as 


y  ( —  1 » T)  =  0 

(8a) 

y  ( 1 ,  r)  =  0 

(8b) 

g(-l.T)  =2X 

(8c) 

gd,T)  =  -2X 

(8d) 
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For  the  case  of  impulsive  loading,  Q  =  0  in  Eq.  (3).  Further, 
taking  account  of  the  initial  shape  of  the  arch,  Eqs.  (3)  and  (4)  may  be  put 
in  the  final  form: 


+ 

dr2  dx4 


92  &L 
dx2 


>2  - 


8 

3 


X2  - 


(9) 


(10) 


To  find  the  static  equilibrium  configurations  we  set  the  inertia 
term  on  the  left-hand  side  in  Eq.  (9)  equal  to  zero  and  get 

iiY  4  Bs  £i=  0  .  ill) 

dx4  dx2 

We  would  like  to  show  nonexistence  of  a  real  number  B2  and  corresponding 
displacement  y(x,8)  ,  other  than  the  trivial  solution  8=0,  y(x,  0)  =  y^  , 
which  satisfy  both  Eqs.  (10)  and  (11)  together  with  boundary  conditions  (8a-d). 

We  first  prove  that  B2  cannot  be  negative.  Multiplying  both  sides  of  Eq.  (11) 
by  y  and  integrating  between  the  limits  from  -1  to  +1  ,  we  obtain 
1  1 

— ^  y  dx  +  3 2  f  — — ^  y  dx  =  0  (12) 

dx4  J j  dx2 


Integrating  the  above  equation  by  parts,  and  making  use  of  boundary 
conditions,  Eqs.  (8)  ,  we  obtain  the  Rayleigh  Quotient: 
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The  right  hand  side  of  this  equation  is  positive  definite  for  any  real¬ 
valued  -  displacement  y.  Hence,  B  can  only  be  real-valued.  The 
physical  meaning  of  02  >  0  is  that  there  cannot  exist  additional  equil¬ 
ibrium  positions  under  tensile  axial  force. 

The  formal  solution  to  Eq.  (11)  can  be  easily  found,  if  one  notes 
(cf.  the  quasi- static  analysis  in  [4]  )  that  in  the  derivation  of  the  shallow 
arch  equations  it  is  assumed  that  the  axial  force  is  constant.  Thus,  if 
8  is  an  arbitrary  real  number,  where  Eq.  (10)  is  considered  as  a 
constraint  on  y(x,  3)  ,  then  the  general  solution  of  Eq.  (11)  can  be 
written  as 

y(x,  8)  =  Aj  sin  8  x  +  As  cos  8  x  +  A3  x  +  A4  .  (13) 


The  unknown  coefficients  A.  ,  i  =  1  ....  4,  are  found  from  the  boundary 
conditions.  Use  of  Eqs.  (8)  results  in  the  following  equations: 


-  sin  8  Ax  +  cos  g  As  -  A3  +  A4  =  0  (15a) 

sin  8  Ax  +  cos  g  A2  +  A3  +  A4  =  0  (15b) 

B  cos  8  +  8  sin  8  A2  +  A3  =  2X  (15c) 

8  cos  8  Ax  -  8  sin  8  A2  +  A3  =  -  2X  .  ( 1 5d) 


Linear  combination  of  these  equations  yields  the  following  values  of 
A2  and  A4  : 


•r»S 


2X 

B  sin  8 


_  2\  cos  8 

4  ~  8  sin  8 

Substitution  of  A2  and  A4  into  Eqs.  (15a)  and  (15c)  gives  the  relations 
for  the  determination  of  Ax  and  A3  as  follows: 
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sin  B  -Ai 


0 


(16a) 


-  A  3  - 

B  cos  B  Aj  +  A3  =  0  .  (16b) 

There  are  only  two  cases  to  be  considered.  If  Eqs.  (15a,b,c,d) 
have  a  unique  solution,  then  the  solution  to  Eqs.  (16a,  b)  is  Ax  =  Aa  =  0  . 
Otherwise,  when  Ai  and  A3  are  nonzero,  the  solution  to  Eqs.  (15a,  b.c,'1) 
is  nonunique,  and  the  vanishing  of  the  determinant  of  Eqs.  (15a,  b,  c,d) 
implies  either  that  B  =  0  ,  in  which  case  the  arch  is  undeformed,  or  that 
the  determinant  of  Eqs.  (16a,  b)  vanish,  viz. 

B  =  tan  0  (17) 


In  this  case 


A3  =  -  sin  B  Ax 


(18) 


CASE  1,  Aj  =  A3  =  0 

In  this  case,  Eq.  (14)  takes  the  form: 

y  =  B  sin  b'  *C°S  PX  "  C°S  '  (19) 

Equation  (19)  is  the  formal  solution  for  the  symmetric  static  equilibrium 
modes,  with  8  to  be  determined  from  constraint  condition  (10).  By  sub¬ 
stitution  from  Eq.  (19)  into  Eq.  (10),  a  transcendental  equation  for  the 
determination  of  B  is  obtained. 

—  B3  ( — —  -  |  }  S  -  4  cot  B  =  0  .  (20) 

X2  \  sinaB  3/ 

The  left-hand  side  at  Eq.  (20)  is  an  even  function  of  B  ■  Hence,  it  is 
sufficient  to  consider  positive  values  of  8  and  show  that  Eq.  (20)  does 
not  have  any  real  roots.  Since  cot  8  is  positive  only  for 
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n  =  0,1,2,... 


(21) 


n  ir  £  8  £  (n  +  j)  U 

we  need  only  consider  roots  in  this  region.  Moreover,  we  will  demonstrate 
that  it  is  necessary  only  to  show  that 

—  g3  +  (— —  -  -|W  >  4  cot  p  0  <g  s  I  .  (22) 

X2  '  sin3  g  3 1  L 

First  note  that  at  g  =  -|  the  inequality  obtained  holds.  Also  for  g  =  ^  , 

3 

cot  g  =  1  and  sin  g  =  1 1  ,  Eq.  (22)  becomes  ll  X2  (-~|)  +  4ir/3  >4  . 

On  the  other  hand  for  small  g  ,  we  can  expand  sin2  g  and  cot  g  in  series 
of  g  and  retain  the  terms  up  to  order  five.  Performing  the  expansion  and 
rearranging  Eq.  (20)  we  get  a  new  inequality  which  is  certainly  satisfied. 

|  —  g6  +  0  (g7)  >0  for  g  sm*ll 

Hence,  inequality  (22)  holds  for  0  <  g  5  tt/2.  Furthermore  it  is 
observed  that  this  inequality  still  holds  if  8  is  replaced  by  g  +  ntT  with 
n  =  1,2,  .  .  .  .  Hence  (22)  holds  for  all  g  in  the  intervals  (21),  and  so 

Eq.  (20)  possesses  no  real  roots. 

CASE  2.  Ai  and  A3  0  ,  g  ~  tan  g 

For  this  case  Eq.  (14)  becomes 

2  > 

Y  =  « — 'A—z  (cos  g.  x  -  cos  g.)  +  Ai  (sin  g.  x  -  x  sin  g.)  .  (23) 

7  g.  sin  g.  1  1  1  1 

Here  the  8  .'s  are  the  eigenvalues  of  Eq.  (11)  determined  by  the  relation 
Bi  =  tan  gj  ,  i  =  1, 2,  3,  .  .  . 


8 


The  constant  Ai  in  Eq.  (23)  is  obtained  through  use  of  the  constraint 
Eq.  (10).  Substituting  Eq.  (23)  into  Eq.  (10),  we  obtain,  after  some 


simplification, 


A,2  = 


» ?  ♦  4  ** 

sin2  B.  -  8? 


For  the  non-existence  of  real  Ai  ,  it  is  sufficient  to  note  that 
the  denominator  to  Eq.  (24)  is  always  negative  for  3.  >  0  . 
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3.  DISCUSSION 


The  conclusion  derived  from  the  foregoing  analysis  is  in¬ 
tuitively  reasonable.  Physically  speaking,  we  want  to  know  if  it  is  possible 
to  deform  the  arch  to  a  certain  shape,  such  that  after  removal  of  the  load 
it  will  not  spring  back.  For  hinged  arches,  there  are  at  least  three  symmetric 
static  equilibrium  states  under  zero  load, and  these  are:  the  original  un¬ 
deformed  shape;  a  nearly  flat  position;  and  a  nearly  inverted  position. 

On  the  other  hand,  if  clamped  arches  are  deformed,  this  deformation  gives 
rise  to  resisting  moments  at  the  ends  which  tend  to  return  the  arch  to  moment- 
free  state.  Thus  it  seems  plausible  that  no  other  static  equilibrium  con¬ 
figuration  exists.  The  same  reasoning  could  also  extend  to  the  case  of 
shallow  shell  structures  which  exhibit  snap- through  phenomenon,  such  as 
clamped  conical  and  spherical  caps.  The  nonexistence  of  snap-through 
under  impulsive  loading  in  shallow  spherical  caps  has  been  discussed  in 
detail  by  Huang  [14]  . 

In  the  light  of  the  present  findings,  it  is  therefore  of  interest 
to  examine  the  works  of  Simitses  and  Humphreys  [1, 2,  7]  ,  in  an  attempt 
to  explain  their  conclusion,  which  is  in  contradiction  to  the  result  of  the 
present  paper. 

Simitses,  in  his  analysis,  represents  the  deflected  shape  of 
the  arch  by  two  symmetric  modes  and  one  antisymmetric  mode  in  the  form 

y  -  yQ  +  (bp  +  bx  cos  £  +  b2  cos  2  £)(1  -  cos  2  £)  (25) 

in  which,  b.'s  are  unknown  constants,  and,  in  terms  of  the  notation  used 
here,  £  =  ^(x  +  1)  and  ^  denotes  the  approximation  to  the  true  deflection 
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y  .  Substitution  from  Eq.  (25)  into  the  expression  for  potential  energy 
and  its  minimization  gives  rise  to  a  set  of  simultaneous  equations  in  terms 
of  unknowns  b.'s.  The  solution  of  these  equations  then  gives  the  static 
equilibrium  points  in  terms  of  generalized  coordinates  b.'s.  It  is  found 
that  there  are  three  symmetric  equilibrium  points  for  this  approximate 
solution,  this  can  be  seen  by  setting  Q  =  in  Eq.  (61)  of  Ref.  [1]  . 

It  is  not  clear  how  Simitses  could  have  found  upper  and  lower  bounds 
for  critical  impulse,  based  on  the  existence  of  symmetric  equilibrium 
positions  alone. 

We  could  cite  two  reasons  as  to  why  Eq.  (25)  could  give  spurious 
results  with  regard  to  determination  of  static  equilibrium  poinls.  First, 
as  is  well  known,  the  use  of  an  approximate  deflection  function  in  con¬ 
junction  with  minimization  of  potential  energy  is  equivalent  to  the  addition 
of  constraints  to  the  system.  Such  constraints  may  be  sufficient  to  give 
rise  to  the  appearance  of  otherwise  non-existent  equilibrium  points.  A 
second,  and  perhaps  more  important  reason  is  the  following.  Referring 
to  Eq.  (14)  and  (21)  ,  it  may  be  noted  that  static  equilibrium  shape  is 
solely  a  function  of  the  axial  compressive  force  in  the  arch  .  This 
axial  force  in  turn  is  given  by  Eq.  (10).  It  can  be  seen  that,  evaluation 
of  g2  is  equivalent  to  finding  the  area  under  the  curve  obtained  by  squaring 
the  slope  of  the  deflection  curve.  Thus,  even  if  Eq  (25)  is  a  reasonable 
approximation  to  the  deflection  shape,  its  derivative  is  less  of  an  approxi¬ 
mation  and  could  introduce  a  large  inaccuracy  m  the  determination  of  the 
axial  force,  thereby  giving  rise  to  the  existence  of  fictitious  equilibrium 
positions. 


11 


It  is  worth  noting  that  the  seconc  effect  mentioned  above  may 
not  be  as  pronounced  in  the  case  of  step  loading.  This  is  because,  in 
such  problem,  the  static  equilibrium  shape  is  a  function  of  external  loading 
as  well,  and  the  relative  inaccuracy  in  B2  could  be  expected  to  be  small. 

Humphreys  has  two  papers  on  the  subject  under  discussion.  In 
[7]  he  appioximates  the  arch  deflection  shape  by  two  symmetric  modes 
of  polynomials.  All  the  objections  applied  to  Simitses  results  are  equally 
applicable  here.  Furthermore,  as  he  points  out,  values  of  the  shallowness 
parameter  must  be  large  enough  for  a  saddle  point  to  exist.  From  the  argu¬ 
ments  4  resented  here,  it  is  clear  that  as  the  shallowness  parameter  becomes 
larger,  the  inadequacy  of  two  mode  representation  becomes  more  apparent. 

In  another  paper  [2]  ,  Humphreys  begins  by  assuming  for  the  de¬ 
flection  shape  a  series  of  the  clamped  beam  eigenfunctions,  from  which 
he  retains  first  six  terms.  Then,  applying  the  Galerkin  technique,  a  set 
of  simultaneous,  nonlinear,  second  order  equations  in  time  is  obtained. 
Solutions  are  obtained  with  an  analog  computer.  For  a  proper  interpretation 
of  the  results,  it  is  necessary  to  discuss  the  procedure  for  determining  the 
critical  snap- through  impulse  (or  load)  in  numerical  calculations. 

The  procedure  for  determination  of  critical  snap- through  impulse 
(load)  is  to  carry  out  deflection- time  solutions  for  different  values  of  load 
parameter.  Then,  by  plotting  curves  of  the  maximum  deflection  (or  maximum 
average  deflection)  against  the  load,  the  point  at  which  a  jump  in  deflection 
occurs  represents  the  critical  value  of  the  load. 

For  the  purpose  of  illustration  of  this  procedure,  let  us  consider 
the  load-deflection  curves  shown  in  Fig.  2  and  3.  These  curves  are  ob- 
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tained  from  a  finite  difference  solution  of  the  equations  for  an  elastic, 
shallow,  clamped,  circular  arch.  The  arch  is  subjected  co  a  suddenly 
applied  cenir<il  concentrated  load  of  constant  magnitude  and  infinite 
duration.  Study  of  these  figures  show  that  the  jump  phenomenon  referred 
to  above  is  quite  distinct  when  the  value  of  rise  parameter  \  is  large. 

For  smaller  values  of  >„  ,  the  jump  disappears  and  is  replaced  by  an  in¬ 
flection  point.  The  appearance  of  an  inflection  point  is  then  an  indication 
that  arch  no  longer  snaps- through. 

In  Humphrey's  analog  computer  solution  [2]  ,  the  impulse-de¬ 
flection  curves  exhibit  only  inflection  points.  In  the  light  of  the  example 
given  above,  existence  of  an  inflection  point  would  not  be  evidence  of 
snapping.  Hence,  the  error  here  is  either  due  to  misinterpretation  of 
numerical  results  or  to  an  imprecise  use  of  the  term  snap- through. 
Humphrey's  experimental  results  are  also  subject  to  similar  objections, 
for  in  [2]  it  was  assumed  that  arch  snaps  when  the  amplitude  of  the 
vibrations  exceeds  the  initial  rise  of  the  arch. 

Up  to  this  point,  we  mainly  have  been  concerned  with  non-existence 
of  snap-through  in  clamped  arches.  However,  indications  are  the  clamped 
shallow  spherical  shells  subjected  to  impulsive  loading  also  do  not  snap- 
through  [14]  .  This  conclusion  is  supported  by  results  of  Thurston  [12]  . 

A  glance  at  his  static  ioad-defiection  curve  shows  that  at  zero-load  there 
exists  no  static  equilibrium  configuration  aside  from  the  undeformed  position. 
Further  evidence  is  obtained  from  the  study  of  other  dynamic  results  [8,  11] 
which  show  only  existence  of  inflection  points  in  continuous  impulse-de¬ 
flection  curves,  and  this  does  not  necessarily  imply  snap- through. 
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